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Abstract 

In this work, we have investigated conduction through an artificial molecule comprising two 
coupled quantum dots. The question addressed is the role of inter-dot coupling on electronic 
transport. We find that the current through the molecule exhibits step-like features as a function 
of the voltage between the leads, where the step size increases as the inter-dot coupling is increased. 
These step-like features disappear with increasing tunneling rate from the leads, but we find that 
in the presence of coupling, this smooth behavior is not observed rather two kinks are seen in the 
current voltage curve. This shows that the resolution of the two levels persists if there is finite 
inter-dot coupling. Furthermore, we also consider the effects of electron-phonon interaction as well 
as dissipation on conduction in this system. Phononic side bands in the differential conductance 
survive for finite inter-dot coupling even for strong lead to molecule coupling. 



A. Introduction 



Ever since the original proposal of a single molecule rectifying diode by A. Aviram and 
M. A. Ratner l| , there has been great interest in electron transport (both charge and spin) 
at the single molecule level. Theoretical as well as experimental work in this direction has 
led to the promising field of molecular electronics jo ] jsj. With an eye on applications it is ex- 
pected that the understanding of quantum electron transport at the molecular scale is a key 
step to realizing molecular electronic devices yll. Experiments on conduction in molecular 

nn 

junctions are becoming more common, |^|5[ and references therein. On the experimental 
front, the most common methods of contacting individual molecules are through scanning 

nn 

tunneling microscope tips and mechanically controlled break junctions |6][8]. Early experi- 
ments focused on the absolute conduction and on trends such as dependence on wire length, 
molecular structure, and temperature. From a theoretical point of view, investigating elec- 
tron transport in an electrically contacted molecule is a challenging problem. In this system, 
the interaction of the electronic degrees of freedom with the vibrational ones of the molecule 
need to be considered. In addition, there can be further complications arising from the 
electron-electron interaction on the molecule as well as effects of the environment surround- 
ing the system. Most of the formal theoretical work on transport in molecular electronics has 



relied on the generalized master equations approach 
function (NEGF) method I?]. 



9| 10| and the non equilibrium Green's 



As mentioned above, an important feature that can affect charge transport in a single 
molecule is the coupling of electrons to quantized molecular vibrational motion, phonons 



For transport through a single level quantum dot molecule a 



ot of work has 



12|(l3|. In this 



already been done taking into account vibrational degrees of freedom fll| 
paper, we consider transport through a single molecule consisting of two coupled quantum 
dots in parallel configuration. Electron transport in double quantum dots has been an 
area of active research, |34i] and references therein. In parallel configuration electron from 
the lead can tunnel through either of the two dots. We address the role of finite coupling 
between the dots on the electronic transport. We also take into account the electron-phonon 
interaction as well as the dissipative effects of the environment. It has been established that 
it is important to take into account electron-phonon interaction in the study of transport in 
single and double dot systems 35|]. Here we show that the inter-dot coupling will significantly 
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affect transport wlien a single electron can occupy either dot in the presence of electron- 
phonon interaction. The difference in the transport properties with and without inter-dot 
coupling will be discussed in detail in this work. For a single level molecular system with 
electron phonon interaction, phonon side band peaks start disappearing with increasing 
tunneling rate|2j] whereas we find that for a coupled dot molecule the phonon peaks survive 
even if the tunneling rate from the leads is increased. 

B. Model 

Our system is a laterally coupled double quantum dot. It is assumed that only a single 
level in each dot participates in transport. We allow finite coupling between the single 
electronic levels of the two dots. An electron from the leads can tunnel through either of 
the two dots. 

The full Hamiltonion describing our system is 

H = Hm + H Leads + Ht- 

It is the sum of the electron Hamiltonian of the coupled dot molecule Hm, the Hamiltonian 
of the leads H Leads, the tunneling Hamiltonian Ht describing the molecule-to-lead coupling. 
We explain each term in the full Hamiltonian separately: 

2 2 

Hm = ^iadl^dio- + '^i,j,o-dif^djf^. (1) 

The first term represents two discrete energy levels, one in each dot. d\^, di„ create and 
annihilate an electron in state \ia > on the dot. The second term represents inter-dot 
coupling where tij^a represents coupling between the electronic states of the two dots. We 
have assumed Uj^a = ij,i,a = t. 

H Leads ~ ^ ^ ^ukaCj^j^^C^ka ■ (2) 
k 

This represents the leads Hamiltonian. Indices u,k,a refer to the left/right leads, the elec- 
tronic wave vector in either lead, and leads electrons spin. 

The tunneling Hamiltonian describes hopping between the leads and the molecule. Direct 
hopping between the two leads is neglected: 

Ht = '^iVi^,y,k,aClkadia + he). (3) 
i,k 
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The first term represents creation of electrons in tlie lead and annihilation of electrons in 
the coupled dots, while the second term represents creation of electrons in the coupled dots 
and annihilation in the lead. Here Vi^^^k,a- denotes lead-system coupling (hopping) amplitude 
and he denotes hermitian conjugation. We consider contacting the coupled dots with two 
metallic leads. 

Finally, phonons, electron-phonon coupling, heat bath and phonons coupling are de- 
scribed by the following Hamiltonian jl?! (In this calculation we work with ^ = 1): 

2 

Hphonon+Bath = ^q(^\0'q + Kia{'^\ + + 

q a=l q 

Y^^p^lh + Y.^'iP^^l + + «^)- (4) 

The electron-phonon interaction is included with in the first Born approximation, which 



is resonable when e 
lem was studied in 



ectron phonons coupling is weak. For a single dot molecule this prob- 



2l| 22|, whereas we consider two dots in the molecule interacting with 



phonons of frequency Uq. The first and the third term represents phononic and heat bath 
energy. Here Wg, are phonon and heat bath energies. a)qaq{b^phis) are phonons creation 
and annihilation operators (heat bath creation and annihilation operators). The second term 
represents electron-phonon interaction and is the coupling strength of this interaction. 
The last term represents phonon and heat bath coupling and Nqp is the coupling strength 
of phonon heat bath interaction. 



C. Method 



Our approach is based on the nonequilibrium Green function technique 



29j|30j, which is 



now a standard technique in mesoscopic physics as well as molecular electronics. We follow 



the formulation pioneered by Meir and Wingreen[3l|], Jauho and co-workers |14| |33l| . The 
case of intermediate and strong electron-phonon coupling at finite tunneling rates is the most 
interesting regime but it is also the most difficult. Only the approaches by Flensberg 24|. 



and Galperin et al 



ITjexist, both starting from the exact solution for the isolated system 



and then switching on tunneling as a perturbation [2]. The current from lead p is given by 
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the well known expression 32 ] 

r(e) 



e 

2^ 



(5) 



Here S^-j^^" represents lesser, retarded and advanced self energies of leads and coupled dot 
molecule. S^j'^^,'" = ^k{^juk9p"^'"'i^)^ii'k} where g^'^'°'{e) is the lesser, retarded and advanced 
Green's function of the leads. |3l[|. 



We employ the wide-band approximation, where the self-energy of the coupled dot 
molecule due to each lead is taken to be energy independent and is given by 



[e - Cfc ± irj) 



^ ^ 0=F = -t2nD{e,)V,,kVlj, = T 



Here D is the constant energy density of the leads. Similiarly the lesser self energy can be 
written as 

5^s<,,,(e) = ^r^V(e) (6) 

k 

where F^^s are the tunneling rates (coupling of leads with the molecule) and /(e) is the 
Fermi-Dirac distribution function. 



Now by ernj 
approximations jl 



oying the current symmetrization and line width proportionality 
we obtain 



de[{GUe)-GtAe)}{f'ie)-ne)}]. 



(7) 



By using the equation of motion technique, we work out the coupled dot molecule Green's 
function and find 



-1 



t 



e - Ci 



e - e,- 



-1 



e - e ,- 



t 



e - e,; 



e - e,- 



Employing Dyson's equation, we can write 



-1 



(8) 
(9) 

(10) 
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The differential conductance in the absence of electron-phonon interaction is 



dr _ 



E 



(11) 



If the electron-phonon interaction is included then our total self energy will become 



Total ^ Leads ~l~ ^el+phoncm- 



and in the Hartree-Fock approximation 



-'el+ph — ^Hartree + ^Fock 



(12) 



where 



and 



ji,Hartree 



^^ji,Focki^) 



/°° df 

/OO J / 



D^,{q,e-e')GUe') + D<{q,e-e')GU, 



(13) 



(14) 



Here Dq represents retarded and lesser free phonon Green's function. At this stage if we 
also include the effects of the heat bath then our zeroth order phonon Green's function will 
be 



1 



1 



(15) 



27 «7 
e - + y e + ^9 + Y 

with 7 = 27r |A^g^p5(e — a;^), |17'], represents dissipation in phonon energy due to its 

contact with the heat bath. 
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i(l + n(e))- 



7 



(16) 



with n{e) the Bose-Einstein distribution function. At zero temperature the imaginary part 
of the electron-phonon self energy is 



ImSji^el+phonon — ^ -^q[^jj(^ ~ ^q) ~ + ^q)] 



(17) 
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where the spectral function Aji(e) = — 2ImGJj(e), whereas 



ReE 



EX2 



Re GUz - Uq) - Re CJz + ujq) + 



2a;„ 



with 



Z7 



(18) 



(19) 



In the presence of electron-phonon interaction and the dissipative effects of the heat bath, 
the differential current is 



dv 27r ^ r^, + rf, 



( (S'lj ( ^) ) ^ + Re E ji ,ei+p^onon) + 



+ ImYiji^el+pho 



(20) 



D. Results And Discussions 



The system under consideration is a molecule comprising two coupled quantum dots. This 
molecule is attached to two metallic leads. The primary focus of this study is the role of inter- 
dot coupling in electronic transport. In addition, the dissipative effects of the environment 
is taken into account by treating the transport in this system surrounded by a heat bath of 
phonons. The vibrational states of the molecule and its impact on electronic transport is 
included through the electron-phonon interaction. Our results show characteristic non-ohmic 
behavior in the current- voltage results presented in Fig.(l). At this stage, we are ignoring 
the phononic and heat bath effects in order to highlight the role of inter-dot coupling. When 
coupling of the system to the leads, which enters through the tunneling-rate, is very small 
then energy levels of the two dots are sharply peaked. Hence, no current flows until the 
apphed voltage bias is in resonance with the level of either of the two dots. We see in Fig.(l) 
that intially no current flows on increasing the bias voltage. But as the applied bias comes in 
resonance with the level of the first dot a sharp increase in current occurs. Further increase 
in applied bias does not lead to increase in current because the level of the second dot is not 
in resonance. As the applied bias is further increased, it comes in resonance with the level of 
the second dot leading to an abrupt increase in current. This explains the step- like features 
with sharp steps and plateaus observed in Fig.(l). Now we consider the effects of inter-dot 
couphng. As a result of the coupling, the levels of the two dots are pushed apart. This 

7 



0.5 1 1.5 2 

Applied Bias (meV) 

FIG. 1: Current as a function of applied bias. The energies of two levels are = O.bmeV, 
= O.SmeF . The tunneling rates from the two states are Ta = 0.04mey, Tij = OmeV, where 

results in the plateaus becoming wider as this requires higher bias voltage before the Fermi 
level in the lead is in resonance with the higher level of the double dot molecule. This is also 
shown in Fig.(l) as the inter-dot coupling is increased. Instead of the inter-dot coupling, if 
we increase the tunneling rate from the lead, broadening of the electronic states in the two 
dots of the molecule takes place. In this situation, current increases linearly with applied 
bias and the system exhibits ohmic behavior. Now if we increase the inter-dot coupling, 
step like features again begin to appear as the applied voltage is tuned since the coupling 
pushes the two levels apart. For sufficiently large tunneling rate the two states broaden to 
the extent that they merge and we find that the current increases smoothly with increasing 
applied bias without any step-like features. Fig. (2). 

These results are also presented in Fig(3) where it is seen that if the lead to the system 
coupling is small then the energy states of the two dots are sharp and this feature appears as 
two peaks in the differential conductance. As the lead to the system coupling is increased, 
the electronic states get broadened. And for sufficiently large tunneling rate (strong lead 
to system coupling) both the states merge and the peaks in the differential conductance 
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Applied Bias (meV) 

FIG. 2: Current as a function of applied bias. The energies of two levels are = O.bmeV, 
= O.SmeV . The tunneling rates from the two states 

disappear. To observe the effects of inter-dot coupling with in the molecule even as the 
tunneling rate from the lead to the molecule (lead to molecule coupling) is increased, we 
show in Fig. (4) that for finite inter-dot coupling, the peaks in the differential conductance 
persist. The electronic states are broadened due to the coupling of the leads and the molecule 
but when inter-dot coupling is taken into account, the difference in energy between the levels 
increases. This compensates the broadening of the levels and allows the two levels to remain 
distinct. This results in peaks corresponding to the two levels appearing in differential 
conductance inspite of broadening of the levels. Fig. (4). 

At the final stage, we consider the effects of the electron-phonon interaction. These are 
shown in Figs (5) and (6). We find that in addition to peaks in the differential conductance 
corresponding to the two levels of the dots there are peaks due to phonons. These phononic 
peaks (side bands) occur as the electrons can exchange energy with the phonons and con- 
tribute to conductance. To focus on the role of inter-dot coupling on phononic peaks, we 
see that as we increase the (lead to molecule coupling) tunneling rate electronic states are 
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FIG. 3: Differential conductance (dJ/dV) as a function of applied bias and tunneling rates. Applied 
bias and tunneling rates are in the units of meV. The energies of two levels are ei = O.BmeV 
,62 = O.SmeV ,t = OmeV. 

broadened to the extent that phononic side bands are not visible. On further increase in 
tunneling rate the two electronic states broaden and merge into each other, Fig. (5). If we 
include inter-dot coupling, not only the electronic states remain distinct but the phononic 
effects are not lost either. This can be seen in Fig. (6) where peaks appear corresponding to 
the two electronic states as well as the phononic side band peaks. Even with an increase 
in tunneling rate from the leads to the molecule, these features persist in the presence of 
inter-dot coupling. 
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Applied Bias (meV) 

FIG. 4: Differential conductance (dJ/dV) as a function of applied bias and tunneling rates. Applied 
bias and tunneling rates are in the units of meV. The energies of two levels are ei = O.BmeV 
,62 = O.SmeV ,t = O.lmeV. 

To conclude, in this work we have focused on the role of inter-dot coupling with in the dot 
molecule on electron transport. We have considered a coupled dot molecule, with inter-dot 
coupling, attached to two leads including electron-phonon interaction and the coupling of the 
molecule with an environment allowing dissipation of the phonons. We find that including 
inter-dot coupling has profound and important role in transport. The step like featues in 
current voltage characteristics of peaks in differential conductance corresponding to the the 
dot energy levels is lost in the absence of inter-dot coupling when strong coupling to the 
leads is considered. Inter-dot coupling allows the two levels to remain distinct with peaks 
appearing in the differential conductance even when broadening of the levels in the dots 
occur for strong lead to molecule coupling. Furthermore, phononic side bands that appear 
in the differential conductance also persist in the presence of finite inter-dot coupling even 
for strong lead to molecule coupling. 
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FIG. 5: Differential conductance (dJ/dV) as a function of applied bias and tunneling rates. Applied 
bias and tunneling rates are in the units of meV. The energies of two levels are ei = 0.5meV 
,62 = O.SmeV ,t = OmeV , oj = O.lmeV ,7 = 0.013mey. 
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FIG. 6: Differential conductance (dJ/dV) as a function of applied bias and tunneling rates. Applied 
bias and tunneling rates are in the units of meV. The energies of two levels are ei = O.bmeV 
,62 = O.SmeV ,t = O.lmeV , u = O.lmeV ,7 = O.OlSmeV. 
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